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ABSTRACT 

The Hubble constant estimated from the combined analysis of the Sunyaev- 
Zel'dovich effect and X-ray observations of galaxy clusters is systematically lower 
than those from other methods by 10-15 percent. We examine the origin of the 
systematic underestimate using an analytic model of the intracluster medium 
(ICM), and compare the prediction with idealistic triaxial models and with clus- 
ters extracted from cosmological hydrodynamic simulations. We identify three 
important sources for the systematic errors; density and temperature inhomo- 
geneities in the ICM, departures from isothermality, and asphericity. In partic- 
ular, the combination of the first two leads to the systematic underestimate of 
the ICM spectroscopic temperature relative to its emission-weighed one. We find 
that these three systematics well reproduce both the observed bias and the in- 
trinsic dispersions of the Hubble constant estimated from the Sunyaev-Zel'dovich 
effect. 

Subject headings: galaxies: clusters: general - X-rays:galaxies masses - cosmol- 
ogy: observation 

1. Introduction 

Galaxy clusters constitute an important cosmological probe, in particular in determin- 
ing the Hubble constant through the combined analysis of the Sunyaev-Zel'dovich effect 
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(SZE) (ISunyaev &: Zerdovichlll972l ) and X-ray observations. Recent high-resolution X-ray 
and radio obs ervations enable one to construct a statistical sample of clusters for the Hq 
measurement. ICarlstrom et al.l (120021 ) compiled the previous results of 38 distar ice determi- 



nation to 26 different g alaxy clusters, and obtained Hq = 60 ± 3 km s ^ Mpc ^ f lReese et al. 



2OO2I : lUzan et al.ll2004l . but see Bonamente et al. 2006). Despite its relatively large indi- 



vidual errors, the mean value of Ho estimated from SZE and X-ray appears systematically 
lower than those esti mated with other met hods: e.g. Hq = 72 ± Skms^^Mpc"^ from the 
distance to Cepheids (IFreedman et al.ll2001l). and Hp = 73±3kms"^ Mpc~^ from the cosmic 
microwave background anisotropy fISpergel et al.l 120071 ). 



Possible systematic er rors in the Hn measurement frona the SZE have been extensively 



studied by several authors (llnagaki. Suginohara &: Sutdll995l:lKobavashi. Sasaki &: Suto 



Yoshikawa. Itoh fc Suto 1998 : Hughes fc Birkinshawl 



1996 



1998 



Birkinshawl I1999I : IWang fc Fang 



2007h : they have addressed a number of physical sources of possible biases including the finite 
extension, dumpiness, asphericity, and non-isothermality of the intracluster medium (ICM). 
Nevertheless they were not able to identify any systematic error that affects the estimate of 
Ho by 10-15 percent. Therefore it has been generally believed that the reliability of Hq from 
the SZE is dominated by the statistics. Given that, the 10-15 percent underestimate bias 
mentioned above, if real, needs to be explained in terms of additional ICM physics beyond 
the simple models used in previous studies. 



Recently, iMazzotta et al.l (120041 ) have pointed out that the spectrosco pic temperature. 



Tspef,, is systematically lower than the emission- weighted temperature, Te„. iKawahara et al. 



( l2007l . hereafter Paper I) investigated the origin of the discrepancy and found that both the 
fluctuation of density and temperature and non-isothermality cause the difference between 
Tgpec and Te„. They also found that the probability density functions (PDF) of both density 
and temperature are well approximated by the log-normal function. 

The aim of the paper is to revisit the origin of the bias based on an observable quan- 
tity, Tspec and the log-normal description for fluctuations of ICM. The paper is organized as 
follows. In §2, we briefly review the conventional method to estimate Hq from the spherical 
isothermal (3 modeling of galaxy clusters. Then we describe several possible sources of the 
systematic bias based on the log-normal description of the fluctuations and the spectroscopic 
temperature. We propose an analytical model for the bias in §3. Non-spherical effects are 
considered in §4 on the basis of triaxial model clusters which include the log-normal fluctu- 
ation and the temperature profile. Section 5 explores the validity of our analytic model for 
the systematic bias using clusters extracted from cosmological hydrodynamic simulations. 
Finally, we summarize our conclusions in §6. Appendix A describes the semi- analytic dis- 
tribution function of the bias of the estimated Hubble constant due to the asphericity of 
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clusters. 



2. Estimating Hq from the SZE in the spherical isothermal l3 model 



A conventional estimate of Ho from the SZE is based on the assumptions that the 
gas temperature is isothermal, T(r) = (= const.), and that the gas density follows the 
spherical /3 model: 

2" 



n{r) = no 



1+ ( - 



1 -3/3/2 



(1) 



where no is the central density, Tc is the core radius, and f3 is the index characterizing 
the density profile. These approximations are insufficient to model the full complexity of 
real galaxy clusters. It has been (implicitly) assumed that the average over a number of 
clusters should significantly reduce the resulting error in the estimate of Ho- While we 
will quantitatively argue below that this is not the case, we su mmarize here the commonly 
adop t ed estimator for Hp in the spher ical isothermal /3 model (llnagaki. Suginohara &: Suto 
19951 : iKobayashi. Sasaki fc Sutdll996l ). 



In this idealistic model, the X-ray surface brightness and y-parameter of SZE at an 
angle 6 from the center of cluster are given by 



A.(Td)ngreG(/3) 
47r(l + 2)4 

npg-T kTci rcG{(3/2) 



1 + 



e 



-3/3+ i 



(2) 
(3) 



where rrie is the electron mass, k is the Boltzmann constant, c is the speed of light, is the 
Thomson cross section, Ax(T) is the cooling function, z is the redshift of the cluster, and we 
define 



G(/3) 



r(3/? - 1/2) 



(4) 



with r(x) being the gamma function. 

Combining equations ([2]) and ([3]), one can eliminate no and estimate the core radius as 

y(0)2 myK{%{) G{(3) 



r c,iso/3(^cl) 



^.(0)47r(aTfcTei)2(l + ^)4[G(/5/2)]2' 



(5) 
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where S'x(O) and y{0) denote the values at 6* = 0, the hne-of-sight through the center of the 
galaxy cluster. Note that the right-hand-side of equation (j5]) is written entirely in terms of 
observable quantities. 

Equation ([5]) corresponds to the estimate of the core radius along the line-of-sight. If 
the cluster is spherically symmetric, it should be equal to the core radius in the plane of the 
sky. With the assumption, the measured angular core radius, 9c fit, is related to the physical 
core radius simply by 



?^c,fit = OcfitdA^z) 



(6) 



with dxiz) being the angular diameter distance of the cluster at z. Equatio ns (151) and (E 
may be combined to estimate the angular diameter distance to the cluster (ISilk fc White 
1978h : 



d 



A,est[Z) 



c,iso/3 



9. 



c,fit 



(7) 



If one obtains ciA,est(^) for a number of clusters at different redshifts, one can estimate 
cosmological parameters by fitting to the angular diameter distance vs. redshift relation, 
dA^z). In what follows, howe ver, we consider the above r aetho dology for the purpose of 
estimating Hq. Thus following llnagaki. Suginohara fc Sutol (119951 ). we introduce the ratio of 
the estimated to the true value of Hq: 

dx -freest fi_ .„^ 

H " V- ^ > 

-no.true ' II 



d 



A, est 



Equations dSD and ^ provide commonly used estimators for the radius of clusters along and 
perpendicular to the line-of-sight, ry, r^, respectively, but they are model-dependent and 
ill-defined for generic non-spherical clusters. We will come back to this issue below (§4 and 
5). Note that /^^ > 1 (< 1) corresponds to over(under)-estimate the true Hq. 

Given the approximations underlying the spherical isothermal (3 model, it is not surpris- 
ing that fn for a individual cluster deviates from unity. A more relevant question is whether 
the average over a number of clusters, (/h), is still systematically larger or smaller than 
unity. If such systematic errors exist, can we correct for them by identifying their physical 
origin ? This is what we address in the present paper. 

In fact there are several prev ious attempts toward the sarne goa l , mainly utilizing numer- 
ically simulated gala. xy clusters (llnagaki. Suginohara fc Sutol ll995l : lYoshikawa. Itoh &: Suto 
19981 : lSulkanenlll999l ). They concluded that departure from the sphericity and the isother- 
mality of clusters results in /^^ 7^ 1, but after averaging over a sample the systematic errors 
are relatively small, | {fu) — 1 1 ^ 5%. Our analysis below is different from the previous ones in 



- 5 - 



adopting the spectroscopic temperature, Tgpec, for Td. Indeed Td is a somewhat ambiguous 
quantify for actual clusters (not isothermal). It has been common in this field to assume 
that the emission-weighted temperature: 



T 



/ n^A^{T)TdV 
Jn^A^{T)dV ' 



(9) 



is approximately equal to Tgpec (the above integration is carried out over the entire cluster 
volume). Thus the previous conclus i on is entire ly based on t he as sumption that Td = Tc„. 
Recently, however, iMazzotta et al.l (120041 ) and iRasia et al.l (120051 ) pointed out that Tspcc, 
estimated by fitting the thermal continuum and the emission lines of the X-ray spectrum, is 
systematically lower than Tew. Furthermore in Paper I we found that the difference between 
Tew and Tgpec could be explained through an analytic model of the temperature profile and 
inhomogeneities in the ICM. We will evaluate applying the model and then comparing 
the numerical simulations in the subsequent sections. 



3. Analytic modeling of systematic errors of Hq for spherical clusters 

Identifying possible systematic errors in the estimate of Hq for realistic clusters is in- 
evitably complicated. In order to address the problem as analytically as possible, we consider 
spherical clusters that follow a density profile of equation ([T]) and a polytropic temperature 
profile but with log-normal density and temperature fiuctuations. While the approach in this 
section is not entirely generic, it is useful in understanding the physical origin of systematic 
errors. The present analytic modeling will be tested against numerically generated triaxial 
cluster samples in §H and against those from cosmological hydrodynamic simulations in ^ 

Our task here is to derive analytic expressions for more general cases, which correspond 
to equations ([2]) to ([5]) in the case of the isothermal j3 model. Let us consider first the effect of 
inhomogeneities in ICM. The X-ray surface brightness at the center of the cluster is written 
as an integral over the line-of-sight: 

^.(0) = ^^^^\^y I n{TfA^{T{r))dr. (10) 

Paper I found that the fluctuation fields defined as 5„ = n(r)/n(r) and 5t = T(r)/T(r) are 
approximately independent and follow the r- independent log-normal PDF, -PLN('^n; c"LN,n) 
and -PLN(^n; crLN,T), where crLN,n and ctln.t denote the standard deviations of the density 
and temperature logarithms. The average of equation (ITOl) over many independent line-of- 
sights can then be computed by integrating over the log-normal PDFs. If we further assume 
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that the coohng function, Ax(T), is dominated by thermal bremsstrahlung (bolometric), 
Abremss(^) oc \/T , we Can rewrite equation ffTOl) as 



= ^^^^^^^,4 J ^(^) Abrcmss(T(r))dr. (11) 

On the contrary, their fluctuations do not affect y{0) because the integrand of the ^/-parameter 
is a linear function of both temperature and density. Thus the inhomogeneity effect is well 
described by the factor: 

Xa = exp(a^N,n-f^LN,T/8)- (12) 

The polytropic temperature profile is expressed as 

T(r) = To [n(r)/no]^-\ (13) 

where Tq is the central temperature (at r = 0), and 7 is the polytropic index. Then we 
obtain 

— '^'^47j-(l + 2)4 J ^^^^ ^hTcmss{T{r))dr 
Abremss(To)ngr,G'(/3(7 + 3)/4) 

47r(l + 2)4 ' ^ ' 

and 

nQa^kTor^G{(3-f/2) 
1/0 = 5 , 15 

respectively. Therefore, the core radius in this model is written as 

1/(0)^ m^c^Ab„(To) G(/3(7 + 3)/4) 
rcpoiyLN ^'^5^(o)47r(aTA;To)2(l + ;2)4 [G(/?7/2)]^ ' ^ ^ 

If one attempts to fit the X-ray surface brightness profile under the assumption of the 
isothermal {3 model, the fitted value of the [3 parameter should be 

& ^ (17) 

since Abremss[7'(r)]n(r)^ oc T{ry/'^n{rY oc [n(r)('^+^)/^]^. In addition, the fitted temperature 
should be equal to the spectroscopic temperature Tgpec. Thus the estimated core radius is 
given by equation ([5]): 

y(0)2 m2c^Ab„(Tspec) G(/3fit) 



c,iso/3 



(-^spec) 



5.(0) 47r(aTA;Tspec)2(l + zY [G{I3^,/2)Y 
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Therefore the systematic bias in the estimate of the Hubble constant in this particular model 
should be 

. _ r,,polyLN _ Abrem..(To)/r2 G{f5{^ + ?,) / A) [G{(3^, / 2)f 

jH,polyLN|iso/3 — 77^ T — Xo 



rc,iso/3(Tspec) " Abrems.(Tspec)/T2p,, /2)Y G(A 



fit J 



•2 



*-bremss 
A-bremss (^spec) /T!,pgj. 



G(/3(7 + 3)A 



2 



X<7XT(Tspec), (19) 



G(/57/2) 

where we define xt that expresses the effect of the temperature structure in the ICM. 
It may be more instructive to rewrite equation f|T9|) as 

r (rp \ XT(Tspec) /'On^ 

jH,polyLN|iso/3 — Xa XTUewJ ,rr \ ' l^^J 

since Td was often assumed to be equal to Tg^- Equation (l20l) makes it clear that the system- 
atic bias in the estimate of Hq results from three major effects; Xo- due to inhomogeneities 
in the ICM, X'viTcw) representing the temperature structure assuming that Td = Tew, and 
finally Xt(Tspoc) / Xt(T'ow) coming from the difference from the spectroscopic and the emission- 
weighted temperatures of the ICM. 

Those three factors can be expressed in an approximate but analytic fashion as follows. 
If we adopt the log-normal PDF for the density and temperature inhomogeneities in the 
ICM, Xa = 6xp(crLN „ ~<^hN t/^) (eq-[I2])- As shown in Paper I, cosmological hydrodynamic 
simulations indicate that ctln,™ ~ 0.2-0.5 and ctln.t ~ 0.2-0.3. Thus Xo- ~ 1.04-1.3. The 
second factor can be estimated by using the analytical relation of Tq and Tew in the current 
model: 

Tew/To = eyip{al^T/2)J{p,'y,rc/r^ir), (21) 

where we assume that the cluster has a finite extension and n{r) =0 for the radius r beyond 
the virial radius of the cluster, r^ir, and we define 

J(3^r/r ) = 2^1(3/2, 3/j[l + 3(7 -l)/4]; 5/2; -r^Jrg) 

with 2Ti(a, /3; 7; () being the hyper-geometric function (see §3 of Paper I). Just for simplicity, 
we neglect the term, exp (cTln t/^), that represents the temperature inhomogeneity because 
it is relatively small for (Jln.t ~ 0.2 — 0.3. If we further adopt Ax(T) = Abrcmss(T) oc -\/T 
then XT(T'ew) reduces to 
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(23) 



The result is plotted in Figured] for typical values of the parameters, and indicates that 
XT(Tew) ranges from 0.8 to 1.0 for 13 = 0.5- 0.8 and 7 = 1.1- 1.2. 



Similarly the third factor can be approximated as 



Xspec- 



spec / 



spec 



A 



brcmss I ew 



XT(Te, 



-^cw -^bremss (T^ 



spec J 



spec 



T 

J. pi 



1.5 



(24) 



Several studies confirmed the systematic underestimate of the spectroscopic temperature 
relative to the emission- weighted temperature, T^^p^./Tp^ = 0.8 — 0.9 from cosmological hy- 
drodynamic simulations (IRasia et al.ll2005l : iKay et al.ll2006l . paper I). If Tgpec/Tew = 0.8 (0.9), 
for instance, Xspec-ew amounts to 0.7 (0.85). 
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Fig. 1. — The bias of Ho due to the temperature profile (eq J23ll as a function of /3. Dashed 
and solid curves correspond to the cases that core radius is 10 percent and 2.5 percent of the 
virial radius, respectively. Red and black colors indicate the polytropic index 7 = 1.2 and 
1.1, respectively. 
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4. Numerical modeling of systematic errors of Hq for inhomogeneous and 

triaxial clusters 

So far we have only considered spherical clusters. Asphericity is definitely another 
important source of error in the estimate of Hq. The errors are expected to be significantly 
reduced by averaging over a statistical sample of clusters randomly oriented with respect 
to our line-of-sight. Nevertheless, if clusters preferentially take either prolate or oblate 
shapes, for instance, the residual errors may not be entirely negligible. This is why we 
addre ss the effect of asphericity on the basis of the triaxial approxirnation for the cluster 



ICM f lHughes fc Birkinshawlll998l : I.Ting fc Sutol 12002 : Ihee fc Sutoll200.i I2004J ). 



To investigate quantitatively the combined effects of gas inhomogeneity and asphericity, 
we numerically create three sets of cluster samples and perform mock observations of Monte- 
Carlo realizations. The first one is spherical, but include random gas density and temperature 
fluctuations according to the log-normal distribution. The second one is triaxial without the 
fluctuations. Both the flrst and the second samples assume isothermality (7 = 1). The third 
one corresponds to the model described in §3 except for added asphericity; the log-normal 
fluctuations, the polytropic temperature structure, and the triaxiality are included. We caH 
them model clusters in order to distinguish them from simulated clusters extracted from 
cosmological hydrodynamic simulations (§5). 

Each model cluster is constructed on (512)'^ grid points within 6 h~^Mpc cubic region 
around the center. We first create spherically symmetric clusters with the gas density profile 
following equation ([1]) with jS = 0.65, = lOOh^^kpc, and uq = lO^^cm^^. The gas is fidu- 
cially isothermal with T = 5 keV, while we also consider the case of polytropic temperature 
profile with Tq = 7 keV and 7 = 1.2. We then add random fiuctuations of gas density and 
temperature according to the r-independent log-normal distributions. The X-ray emissivity 
is computed with SPEX version 2.0 assuming collisional ionization equilibrium, the energy 
range of 0.5 — 10.0 keV and a constant metallicity O.SZq. Triaxial model clusters are con- 
structed simply by stretching spherical clusters along the three axis directions by a factor of 
Aq, Aft, and Ac, respectively. 

In mock observations, we extract the quantities necessary to compute rc^isopiTd) and 
refit via equations ([5]) and ([6]) in the following manner. We first fit the projected profiles of 
Sx{0) with a functional form S'x(0)[l -|- (r/rc^fit,sx)^]~^'^^''''^''^ from r = to r = lh~-^Mpc over 
1024 random LOSs toward each cluster. For each LOS, we also compute y{0) and, unless 
otherwise stated, use it directly in our analysis. We will discuss other choices of obtaining 
y{0) in §5.3. As will be described later, the gas temperature Td is obtained by either fitting 
the mock X-ray spectra or simply using the input temperature, depending on the purpose 
of the analysis. We use the template of the spectral energy distribution computed using 
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SPEX 2.0 assuming coUisional ionization equilibrium, the energy range of 0.5 — lO.OkeV and 
a constant metallicity 0.3^0. Assuming that = r c,iso/3 

(Tci) and = r^fit.Sx, we calculate 

fu for each LOS. 

To quantify the bias due to the projection effect, we also compute the volume-averaged 
radial profile of the gas density, directly from the grid data within the radius lh~^Mpc. 
By fitting the profile to the /? model, we obtain the estimated core radius rc,fit,3D, which is 
independent of LOS. We will compare the values of using r^^ = rc,fit,3D and rc,fit,sx in what 
follows. 
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Fig. 2. — The average and rms of fn of the model clusters; (a) spherical clusters with gas 
inhomogeneities and no temperature gradient, (b) ellipsoidal and isothermal clusters, and (c) 
ellipsoidal clusters with temperature gradient and gas inhomogeneities. Crosses and pluses 
denote fa adopting = Tc,fit,sx and = rc,fit,3D, respectively. Thick horizontal lines indicate 
analytical estimations for Xa (dashed), XaXTiTe^) (dotted), and Xo-XT(T'spec) (solid). In panel 
(c), black symbols indicate using rc,iso/3(T'ci = T^pec), red symbols rc,iso/3(^ci = T^^), and 
blue symbols rc,poiy/3, which correspond to the isothermal fit with T = Tgpec, the isothermal 
fit with T = Tew, and the polytropic fit, respectively (see the main text for details). 



Figure [2|i shows the mean and rms values of fn for spherical clusters with no temperature 
gradient. We consider two cases for the log-normal density and temperature fluctuations 
with (ciLN.n, o"ln,t) = (0.5, 0.0) and (0.5, 0.3). The latter set corresponds to the typical value 
for the simulated cluster as Paper I reported. To present the bias produced solely by gas 
inhomogeneities, we here adopt for the volume averaged temperature instead of fitting 
the mock X-ray spectra. It is evident that the fluctuations in gas density yield fn ~ 1-3 
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(i.e. overestimating of Ho by ~ 30%), while those in gas temperature do not contribute 
significantly to the bias. The mean value of is in good agreement with our analytical 
expectation Xa (dashed horizontal lines). The bias due to the projection is ~ 10%. 

The bias produced by ellipsoidal shapes is displayed in Figure [2}d for the two sets of 
the axis ratio (Aa : : Ac = 0.6 : 0.7 : 1 and 0.9 : 0.9 : 1). These sets are the typical 
value of the simulated cluster. Again, to present the bias solely from asphericity, the gas is 
assumed to be isothermal without any fluctuations and we adopt = Tq in computing fjj. 
The average bias is relatively small (< 15%) and that due to the projection is ~ 3%. These 
results indicate that the bias due to asphericity, after averaging over a statistical sample of 
clusters, is smaller than that from gas inhomogeneity. 

Figure ^ illustrates the bias in a more realistic case; we create ellipsoidal clusters with 
the polytropic temperature profile and fluctuations. Two sets of axis ratio (A^ : A;, : Ac = 
0.6 : 0.7 : 1 and 0.9 : 0.9 : 1) are chosen adopting (o"LN,n, c"ln,t) = (0.5,0.3) and (0.3,0.2), 
respectively. In this panel, we show the values of fn based on the following three methods, 
so as to understand clearly the physical origin of the overall bias in the Hq estimation. 

The first method (black symbols in Fig. ^) corresponds to the most conventional case 
of using the isothermal /3-model and the spectroscopic temperature Tgpec. To obtain Tspec 
we fit the mock X-ray spectra from the central (r < 1 h~^Mpc) region of each cluster using 
XSPEC assuming a single temperature MEKAL model. We assume the perfect response, 
and ignore its effect on the spectral temperature (Paper I). Clearly, the value of Hq is 
underestimated by ~ 10-20%. This is in good agreement with our analytical estimation for 
/H,poiyLN|iso/3 = Xo-Xt(^spcc) for a Spherical cluster (solid horizontal lines). To obtain XT(^spec), 
we compute the volume-averaged profile of the density and the temperature. These profile 
are fitted to equation ([T]) and equation (IT^ taking uq, rc,fit,3D, P, Tq and 7 as free parameters. 
We use the adopted values (o"LN,n5 o"ln,t) = (0.5,0.3) and (0.3,0.2) to compute Xa- 

The secon d method (red symbols in Fig. [2t:) aims to mimic previous nume rical studies 
of the Hq bias (llnagaki. Suginohara &: Sutdll995l : lYoshikawa. Itoh &: Sutdll998l ) and adopts 
the isothermal /?-model and the emission-weighted temperature Te„. We obtain by 
directly summing up the temperature of each grid point from the central (r < 1 h~^Mpc) 
region. Also plotted for comparison is an analytical estimate Xo-XT(T'ew) (dotted lines). 
The values of x<j and XT(^ew) are computed as described above. In this case, Xo- and 
Yt(?"pw) pra ctically cancel each other, and i s clos e to unity, consistent wit h the previous 



findings of llnagaki. Suginohara fc Sutd (119951 ) and lYoshikawa. Itoh fc Sutd (119981 ) . This 



shows that the absence of the bias in previous studies is simply an artifact of using Te„, 
which is systematically larger than Tgpec. 
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The third method (blue symbols in Fig. [2t) attempts to eliminate the bi as due to the 



temp erature gradient by using the polytropic profile to estimate the core radius (lAmeglio et al. 

y{0)' myK{%) G(/37/4 + 3/3/4) 
rc,poiy/3 5.(o)47r(aTfcTo)2(l + ^)4 [G{[5^/2)f ' ^ ' 

where we adopt Tq and 7 from fitting the volume-averaged temperature profile of the model 
clusters. The value of (3 is obtained from equation f|T7|) using (3f^^ and 7. The value of so 
obtained should represent the bias arising from sources other than the spectral fitting and 
the temperature gradient. Given the good agreement with the analytical estimate for x<7 
(dashed lines), we conclude that the bias in this case is dominated by the effect from gas 
inhomogeneities. 

In summary, there are three major sources for the bias of i^o; the spectral fitting, the 
temperature gradient, and local density fiuctuations. The former two leads to an underes- 
timate while the latter an overestimate of Hq. In every case studied here, the bias due to 
asphericity is much smaller than the other three. 



5. Comparison with clusters from cosmological hydrodynamic simulations 

5.1. Cosmological hydrodynamic simulations 

We now compare the bias described in the previous section with simulated clusters. 
They are extracted fro ni the Smooth i ng Pa rticle Hydrodynamic (SPH) simulation of the local 
universe performed by iDolag et al.l (120051 ) assuming ACDM universe with f2om = 0.3, = 



0.04, (Tg = 0.9, and h = 0.7. The numbers of dark matter and SPH particles are ~ 20 
million each within a high- resolution sphere of radius ~ llOMpc, which is embedded in a 
periodic box ~ 343 Mpc on a side that is filled with nearly 7 million low-resolution dark 
matter particles. The simulation is designed to reproduce the matter distribution of the local 
universe adopting the initial conditions based on the IRAS galaxy distribution, smoothed 
over a scale of 4.9/i~^Mpc. We choose the six massive clusters identified as Coma, Perseus, 
Virgo, Centaurus, A3627, and Hydra. Figure [3] shows projected surface density maps of these 
simulated clusters. The values of P and 7 of these clusters are listed in Table [H The cubic 
region of 6 h^^ Mpc around the center of each cluster is extracted and divided into 512'^ cells. 
The density and temperature of each mesh point are calculated from SPH particles using 
the B-spline smoothing kernel. A detailed description of this procedure is given in Paper I. 

We perform mock observations over 1024 LOSs for each simulated cluster in a similar 
manner to §1] except for the following points. First, we compute Tgpec and within the 
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virial radius instead of 1 h ^Mpc. Second, we use the fitted value of o"LN,n and o"ln,t in 
calculating Xu of the analytical model. 

Table 1: Properties of the Six Simulated Clusters. 



/5fit,3D * 7 * Aa/Ac Afe/Ac (rc,fit,Sx)/?^c,fit,3D 



Coma 


0.74 


1.17 


0.59 


0.64 


1.03 ± 0.14 


Perseus 


0.64 


1.09 


0.49 


0.61 


1.04 ± 0.18 


Virgo 


0.60 


1.15 


0.44 


0.61 


1.16 ± 0.31 


Centaur us 


0.69 


1.17 


0.68 


0.78 


1.03 ± 0.13 


A3627 


0.69 


1.15 


0.79 


0.83 


1.08 ± 0.06 


Hydra 


0.70 


1.22 


0.84 


0.93 


1.03 ± 0.05 



* 'i'he values of /:^fit,3D and 7 are slightly changed from that listed in Paper I 
due to the improvement of the routine of fits. 



5.2. Results 

Figure H] displays a set of histograms of fn for the simulated Coma cluster. The same 
analysis is done for the other five clusters. Histograms in different colors correspond to the 
symbols of the same color in Figure [2t, and indeed show similar trends for each component of 
the bias. Since the physical length of clusters along the LOS is not symmetrically distributed 
around its mean, the corresponding histograms of are skewed positively. In Appendix A, 
we compute the distribution for the two extreme cases, the prolate and the oblate ellipsoids, 
and find that they yield positively and negatively skewed distributions, respectively. Indeed 
this is consistent with the fact that the simulated Coma is nearly prolate (Table [1]) . 

The (simple arithmetic) mean, (/^) is plotted in Figure O for six simulated clusters. The 
quoted error bars indicate la standard deviation from the mean. Except for the simulated 
Virgo cluster, (/h) is below unity, i.e., Hq is underestimated. It is remarkable that a simple 
analytical model for systematic effects (solid, dotted and dashed horizontal lines) described 
in §3] can reproduce the bias in the simulated clusters. 

We have made sure that the bias from other sources is minor; first, if a cluster has a 
finite extension and is bounded within the virial radius, the value of {fn) becomes smaller 
by ^ 5% (open circles in Fig. [5]). Second, we compute the axis rat i o (An < A;, < Ac) of 



each simulated cluster, basically following the method of IJing fc Sutd (120021 ). but using the 



gas density not the dark matter density. The isodensity surfaces corresponding to the gas 
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Perseus 

'■ w 




A3627 

«. 




Virgo 


f 


Hydra 
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Fig. 3. — Projected surface density maps of the six simulated clusters. Five different {ue = 
3 X 10~^, 1 X 10"^, 5 X 10~^, 3 X 10~^, and 1 x 10~^ [cm~^]) isodensity surfaces are indicated 
with different colors (red, orange, yellow, green, and blue, respectively). The left panels of 
each cluster indicate the view from our galaxy. The right panels are the projection of each 
simulated cluster as seen by a distant observer located to the "right" of each panel on the 
left. The horizontal yellow lines indicate the physical size of lh~^Mpc. 
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Fig. 4. — Distribution of over 1024 LOSs for the simulated Coma cluster. Black, red and 
blue histograms indicate the results for the isothermal fit with T = Tspec, the isothermal fit 
with T = Tew, and the polytropic fit, respectively. 
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/h ,polyIiN|iso|9 ^(j^T^^spec^ 
X(j-Xx(^ew) ^CT 



Centaurus 



A3627 



Hydra 



Fig. 5. — The average and rms of for the six simulated clusters. Black, red and blue sym- 
bols with error-bars indicate the results for the isothermal fit with T = Tgpec, the isothermal 
fit with T = Tew, and the polytropic fit, respectively. Crosses and pluses denote f„ adopting 
r± = Tc,fit,Sx and = Tcfit.sD, respectively. Open circles indicate fa adopting = Tcfit.Sx but 
assuming that the clusters are not extended beyond the virial radius. Thick green horizontal 
lines indicate analytical estimations for Xa (dashed), Xo-XtI^cw) (dotted), and Xo-XT(^spec) 
(solid). 
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densities of rig = 3 x 10"^, 1 x 10"^, 5 x 10""^, 3 x lO"'^, and 1 x 10""^ [cm"^] are shown in Figure 
[3l After eliminating substructures, the axis ratio is calculated by diagonalizing the inertial 
tensor of each surface. The averaged axis ratios of the five different density regions (Table 
[T]) are similar to the typical value adopted in §4. Therefore, we conclude that the spherical 
approximation itself is not a major source of the bias for the simulated cluster. Finally, the 
bias due to the projection is small (crosses and pluses in Fig. |5]). We list in Table [1] the 
average values of rc,fit,sx over 1024 LOSs relative to rc,fit,3D- The ratio is unity within 10 % 
(except for Virgo that has a relatively large dispersion), and basically all consistent with 
unity within the uncertainty. 

It is interesting to emphasize here that the shape of the distribution of reflects the 
shape of clusters fr om the perspective of measuring the three-dimensional shape of clusters 



(iSereno et al.l 120061 . e.g.). If all clusters had the same shape, the observation of one cluster 
toward multiple directions might correspond to that of multiple clusters toward each LOS. 
Of course, the real shapes vary from cluster to cluster. However, if clusters tend to be prolate 
(oblate) preferentially, the distribution of i?o,est should be skewed positively (negatively) as 
shown in appendix A. Therefore, independently of the knowledge of real value of i?o,truej 
the statistical information about the shape distribution may be obtained in principle by the 
distribution of ifo.est- 



5.3. Comparison with previous studies 



The above results are consistent with the previous results of f„^\ with Td = 
Tp„ (Inagaki. Suginohara fc Suto 1995 : Yoshikawa. Itoh fc SutolE998l ). On the other hand, 
Ameglio et al.l (120061 ) explored the bias of c?a using the cosmological hydrodynamic simula- 
tions, and reported that Hq is overestimated by more than a factor of two if one adopts the 
isothermal (3 model. This is opposite to our conclusion here, and we found that this should 
be ascribed to the sensitivity of /h on the adopted values of y(0), i.e., /h oc (i^est ^ 2/^^(0) 
as explained below. 



Ameglio et al.l (120061 ) obtained ?/(0) by fitting the noise-less profile of y{6) up to R[ 



500 

fixing other /3-model parameters from the X-ray profile, while we use directly the projected 
values of ?/(0) in the sim ulation data. The di fference in these two methods is apparent in 
Figure 5 (right panel) of lAmeglio et al.l (120061 ): their fit (solid line) was affected largely by 
the data points at large radii and yielded a value of ?/(0) smaller by ~ 50% than the actual 
data. This enhances the value of /h by more than a factor of two, and indeed accounts 
for their apparently opposite conclusions. We have checked that other differences between 
their analysis and ours (the use of mass-weighted temperature for and the removal of the 
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cluster central region) do not affect the results significantly. 

As far as the bias in the previous SZE observations is concerned, we believe that our 
method is more relevant to what has been done with the real data, because these observations 
were not capa ble of constraining the radial profile of the y-parameter up to large radii with 



high S/N (see lKomatsu et al.ll200ll : iKitayama et al.ll2004j . for the currently highest angular- 



resolution observation of the SZE). 

We have further made sure that the effects of the finite spatial resolution of the ob- 
servations and the central cooling region, which were neglected in our preceding analysis, 
are minor; first, we have evaluated y{0) by fitting y{6) within a radius of 100 /i~^kpc and 
200 h'^kpc. These approximately correspond to the typical angular resolution of the SZE 
observation (~ 1 minute) at z = 0.1 and z = 0.3, respectively. The values of j3 and 9c are 
fixed from the X-ray profile. For the six simulated clusters, the resulting values of /h differ 
from our initial analysis (see Fig. [5]) by —6% to +7% {+2% on average) for r < 100/i~^kpc, 
and by —3% to -|-10% (+5% on average) for r < 200/i^^kpc. 

Second, we have also performed the fitting separately for the X-ray and SZE profiles. 
The values of ^x, /3 and 9c are evaluated by fitting Sx{9) with equation ([2]), while that of 
y{0) is obtained by fitting y{9) with equation ([3]) independently of the X-ray profile. As a 
result, the values of /h differ from those of Figure [5]by —6% and +1% (—3% on average). 

Note that observationally there are several different ways to evaluate y{0), S'x(O), /5 and 
9c- A conventional method is to fit the X-ray imaging data 5'x(^) first. Then the SZ image 
is fitted to obtain y{0) assuming the values of /? and 9c from the X-ray data. Our an alysis 



procedure adopted here follows the conventional method. While iReese et al.l (120021 ) have 
determined from the joint fit to the X-ray and SZE imaging data, the result is almost 
equivalent to the conventional method since the X-ray imaging data have a much higher S/N 
than the SZE data. 



6. Conclusions 

We considered various possible systematic errors of Hq from the combined analysis 
of the Sunyaev-Zel'dovich effect and X-ray observations. In particular we addressed the 
validity and limitation of the spherical isothermal (3 model in estimating Hq, which has been 
used widely as a reasonable approximation after averaging over a number of clusters. We 
introduced the ratio of the estimated to the true Hubble constant, fn, to characterize the 
systematic errors. We constructed an analytic model for J'h, and identified three important 
sources for the systematic errors; density and temperature inhomogeneities in the ICM, the 
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temperature profile, and departures from spliericity. Except for the non-splierical effect, tlie 
most important analytical expression tliat summarizes our conclusion is equation fl20l) . or 
equivalently, 



Honest 

Ho. 



(26) 



true 



In our analytic model discussed in §3, the inhomogeneity bias, Xcr, the non-isothermality 
bias, XTiTevj), and the temperature bias Xspec-ew are given by equations (fT2l) . (123|) . and fl24|) . 
respectively. 

While the above model prediction is fairly general, the net value of /h sensitively de- 
pends on the degree of the inhomogeneity and multi-phase temperature structure of real 
ICM. Our simulated cluster sample implies that Xo- ~ (1-1 — 1-3), XT(^ew) ~ (0-8 — 1), 



Ysper- ew ~ (0.8 — 0.9), aud therefore (/h) ~ (0.8 — 0.9). Given the result of iReese et al. 



( 120021 ). this is certainly indicative, but may need to be interpreted with caution because the 
result is critically dependent on the reliability of the adopted numerically simulated clusters 
as representative samples of clusters observed in the real universe. Exactly for this reason, 
we are attempting more direct (not statistical) comparison of our model prediction against 
observed cluster samples, which will be presented elsewhere hopefully in the near future 
(Reese et al. in preparation). 
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A. Distribution of for prolate and oblate ellipsoids 

In this Appendix, we derive the distribution of fn due to the asphericity of clusters, 
by considering the following two extreme cases; the prolate (Aa = Af, < Ac) and the oblate 
(Aa < Xb = Ac) ellipsoids. We choose z-axis as the long (short) axis and x- and y-axes as 
the short (long) axes for an prolate (oblate) ellipsoid. The direction of the unit vector along 
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Fig. 6. — Schematic representation of a prolate cluster with an axis ratio : Aj, : Ac = 1 : 
A : A. Given the symmetry around the z-axis (the long axis), an LOS through the cluster 
center is specified by an angle 6'a from the z-axis. An arbitrary position in the cluster r is 
expressed in terms of L (the projection of r onto the LOS direction), R (the projection of r 
onto the plane normal to the LOS), and C (the azimuthal angle on the plane normal to the 
LOS). 
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the LOS of an observer, a, is defined in terms of the spherical coordinate {9a., (pa)- Figure [6] 
shows a schematic picture of a prolate ellipsoid. 

Let us define the quantity A = Xa/Xc = Xh/Xc (A = Xb/Xa = Xc/Xa) for the prolate 
(oblate) ellipsoid. We assume that the gas density follows the prolate (oblate) (3 model: 

n{v)\e^ = no(l + (r7re)2)-'^/% (Al) 
f = Irlisin^^r/A' + cos'^r]'/', (A2) 

and 9j. is the angle between z-axis and r. Because the surface brightness profile is independent 
of (pa due to the 2;-axial symmetry, one can express fn as function of 6a- 

For an isothermal cluster, the surface brightness averaged over a circle of radius R is 
proportional to J n^dL averaged over the angle C, of the circle. We put C = where r is 
located on the same plane defined by the LOS and z-axis. To compute the averaged surface 
brightness, we need an expression of the density n as a function of L, i?, and C,. In the 
Cartesian coordinate, r(^ = 0) is 



yE? + I? cos 0a sin (arctan R/ L + 9a 
r(C = 0)=| + L"^ sin 0^ sin (arctan R/L + Oa) |. (A3) 

\/R^ + LP' cos (arctan R/L + 6a_ 

Multiplying the rotation matrix around a, Ma(C) to r{( = 0), we obtain 

r(C) = Ma(C)r(C = 0) 

R cos ( cos 0a COS 9a + R siu ( siu (j)a + L COS 0a siu 9a 

R COS C sin 0a COS 9a — R sin C, cos 0a + sin 0a sin 9a I , (A4) 

L COS 9a — R COS C sin 9a 

where 

Ma(C) = 

xi+{l-xl)cOsC XaVai^ - COSC) + ZaSmC, ZaXa{l - COS Q - Va^in C, 

XaVai"^ - COS C) - slu ( yl+{l- yl) COS C VaZai^ - COS C) + Xa slu ( | , (A5) 

2aa;a(l - COSC) + 2/aSinC ?/a2;a(l - COSC) - XaSinC zl + {I - Z^) CO?. C, 



and 



Xa\ 




Va 










Thus, we obtain 



(A6) 



|r(^) I cos 9j. = L cos 9a — R cos ( sin 9a- (A7) 
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Then, |r(C)| and Or are written as 

|r(C)| = + i?2 (A8) 

/Lcos6'a - i?cosCsin6'a\ 
c^r = arccos , . (A9) 

Combining with equation (lAip . we can write n(r)|0^ in terms of and C as 



where 



n(r)|9. = n„ (1 + (f (fl, L, Ois./i-c)') = "(B. i, C), (AlO) 



Then, the averaged surface brightness at R is 
^x(^)l^a = dL[niR,L,C)\eJ 



27r />oo 



/ dqL 



oo 



?! + 4 + (^^ - 1) COS - g/j COS C sin 



+ 1 



-3/3 



A2 

- . HQRh. (A12) 

where we define the normahzed length by Vc, qn = Rjrc^qL = L/vc- We compute I{qR)\e^ 
numerically for A = 0.5 (prolate) and A = 2.0 (oblate) adopting /? = 0.65. We fit I{qR)\e^ 
from qji = to qn = 10.0 by with a functional form of the surface brightness profile assuming 
the spherical beta model (oc [1 + {qR/qc,fit\e^)'^]~^^^^^'''^~^^^'^)- Thus, we obtain the counter part 
of '^c.iso/j , '?c,fit|6»a = ''^c,fit,Sx/''^c and the fitted value of (3, (3fit\e^- While, gc,iso/3|ea = ^cAso/s/fc is 
written as 

gc,iso,,a = (sm e./A + cos e.) ^^^^^y^^^^^^^j (A13) 

The first term of the right-hand side represents the elongation of the radius toward the LOS. 
The second term is the correction to the use of Pfit\d^ in observation instead of the true /3. 
However, the correction is very small (within 0.01% error). 

Finally, we obtain the bias of Hq as a function of 6'a, 

UO.) = (AM) 

Q'c,iso/3|6la 

The probability of for the random assignment is proportional to the solid angle Qlfu). 
If /H(6'a) is a monotonic function, the PDF of is obtained as 



de^ 



dfn 



- 23 - 



dfu 



(A15) 



where O^Uh) = f^H 



Dotted lines in the upper panel of Figure [7] show equation (lAlSP for prolate (A = 0.5) 
and oblate (A = 2.0) ellipsoids. As shown in the lower panel, the corresponding 6a is a 
monotonically increasing (decreasing) function of for the prolate (oblate) ellipsoid. At 
= 0, fn is equal to A, which corresponds to the case that the LOS is along the 2;-axis. 



The PDF diverges at 6^ 
flAlSp imply that 

P{U)cxsm9a{U] 



-- 7t/2. This can be understood as follows. Equations flA13P to 



^c,iso/3|^a 



-1 ^ vW0a(/«) + A~2sin^0a(/. 



COsOaifn 



(A16) 



where we ignore the 6'a-dependence of gc,fit|ea P&tie^- Thus 9a ~ 7r/2, Pifn) diverges 
as 1/ cos6'a. Note, however, its integration over a finite size of fn does not diverge (see 
eq. |A15j ). This is plotted in the sohd histograms, where the bin size A/^ = 0.05 is adopted. 
The resulting distribution is skewed positively (negatively) for the prolate (oblate) ellipsoid, 
which is consistent with the results shown in Figure |H 
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